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ABSTRACT 


Some  empirical  relations  between  diffusion  characteristics  are 
investigated  by  the  use  of  carefully  examined  data  from  dye  release 
experiments  in  the  surface  layer  of  the  sea.  These  data  cover  a 
time  scale  of  diffusion  ranging  from  1  hour  to  1  month  and  a  length 
scale  from  100  m  to  100  km.  Two  "oceanic  diffusion  diagrams"  are 
prepared;  one  showing  horizontal  variance  versus  diffusion  time  and 
the  other  showing  apparent  diffusivity  versus  the  scale  of  diffusion. 
The  overall  behaviors  of  the  horizontal  variance  and  of  the  apparent 
diffusivity  are  evidently  different  from  those  which  the  similarity 
theory  of  turbulence  deduces.  However,  there  still  remains  a 
possibility  that  the  similarity  theory  may  be  valid  locally  in  time- 
or  length-scale  with  a  variable  parameter,  i.  e.  ,  the  rate  of 
turbulent  energy  transfer.  The  diagrams  provide  a  practical  means 
to  predict  the  rate  of  horizontal  spread  of  substance  introduced  from 
an  instantaneous  point-source  as  well  as  the  apparent  diffusivity  as 
a  function  of  the  size  of  diffusing  patch. 
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1.  Introduction,  Turbulent  diffusion  of  contaminants  iB  not  only  a 
matter  of  theoretical  interest  but  is  also  a  problem  of  practical 
importance  in  connection  with  the  discharge  of  industrial  wastes 
including  radioactive  materials  into  the  oceans.  The  process  of 
oceanic  diffusion  is  quite  complex,  and  hardly  any  single  theory  can 
explain  or  interpret  the  entire  pattern  of  the  diffusion.  In  addition  to 
theoretical  studies,  therefore,  progress  in  oceanic  diffusion  is  still 
very  much  dependent  on  an  empirical  approach  by  means  of  diffusion 
expe  riments . 

Among  various  types  of  tracer  used  in  experimental  studies  of 
oceanic  diffusion,  fluorescent  dyes  have  proved  the  most  promising, 
especially  Rhodamine  B.  It  has  been  regarded  as  one  of  the  most 
accessible,  stable,  harmless,  and  convenient  fluorescent  tracers 
for  use  in  sea  experiments  (Carpenter,  I960,  Okubo  et  al  1957). 

Since  Pritchard  and  Carpenter  (I960)  developed  a  new  field 
technique  for  the  direct  continuous  observation  of  the  fluorescent 
concentration  with  the  aid  of  a  fluorometer,  experimental  data  of 
oceanic  diffusion  have  been  accumulated.  In  particular,  an  internation¬ 
ally  cooperative  experiment,  operation  "RHENO”,  in  the  North  Sea 
made  it  possible  to  cover  a  scale  of  diffusion  ranging  up  to  100  km 
and  a  time  scale  up  to  one  month.  Toward  the  other  end  of  the  scale, 
a  set  of  experiments  off  the  coast  of  Southern  California  provided  the 
data  of  small  scale  diffusion  of  30  m  or  so.  A  fairly  large  amount  of 
data  have  been  obtained  for  intermediate  scales  of  diffusion  ranging, 
say  from  1  km  to  10  km.  So  far  no  appropriate  d.-.ca  are  available 
from  dye  experiments  in  the  deep  layer.  Hence  we  are  restricted  to 
data  obtained  in  the  surface  layer. 

Using  the  available  data  of  horizontal  diffusion  until  1962, 

Okubo  (1962)  presented  a  set  of  diagrams  showing  the  empirical 
relations  between  the  variance  of  the  concentration  distribution  and 
diffusion  time  and  between  the  coefficient  of  eddy  diffusion  and  the 
scale  of  diffusion. 
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In  this  report,  a  new  set  of  diffusion  diagrams  are  constructed 
on  the  basis  of  new  data  accumulated  since  1961.  None  of  the  data 
chosen  here  was  used  in  the  diffusion  diagrams  constructed  in  1962 
in  order  to  obtain  completely  new  diagrams.  Carefully  examined 
data  were  analysed  by  a  consistent  method  in  order  to  compute  the 
variance  of  the  horizontal  distribution  of  concentration.  One  of  the 
purposes  of  the  present  study  is  to  compare  the  findings  in  the  old 
diffusion  diagrams  with  the  new  set  of  diagrams. 

2.  Data  used  in  the  present  study.  After  exhaustive  search  the 
author  has  found  20  sets  of  diffusion  experiments,  which  provide  data 
suitable  for  the  present  purpose.  The  criteria  for  defining  suitable 
data  are  as  followc  : 

1)  The  release  of  dye  should  be  as  close  tc  a  type  of 
instantaneous,  point  source  as  possible;  the  duration  of  release  and 
the  initial  (horizontal)  size  of  dye  patch  must  be  reported  or  somehow 
estimated  at  least  to  an  order  of  magnitude.  For  a  known  initial  size 
o  -  more  precisely  a  known  initial  variance,  a  critical  time  after  which 
the  diffusion  is  regarded  practically  as  from  a  point  source  can  be 
est^-nated  from  a  theoretical  relationship  for  the  rate  of  change  of 
variance,  e.  g.  ,  the  theory  of  Joseph  and  Sendner(  1958)  gives 

o  r  =  6  P2t  ?  ,  (1) 

o  o 

o 

where  a  :  the  variance  associated  with  the  initial  distribution  of 

o 

patch,  P;  a  diffusion  velocity  after  Joseph  and  Sendner,  and  t^;  a 

characteristic  time  of  diffusion  of  patch  associated  with  the  initial 

size,  o0  .  Obviously  the  critical  time  must  be  much  larger  than  tQ. 

An  arbitrary  choice  of  the  critical  time  herein  made  is  10  times  of 

t  .  Diffusion  data  prior  to  the  cri  cal  time  are  either  disregarded 
o 

or  corrected  to  an  appropriate  diffusion  time  by  adding  the  critical 
time  to  the  actual  diffusion  time. 
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2)  The  dye  patch  should  stay  a  sufficient  distance  from  vertical 
boundaries  so  that  the  field  of  diffusion  may  be  regarded  as  extending 
to  infinity  in  the  horizontal  direction.  In  most  cases  whether  a  dye 
patch  has  felt  the  boundaries  or  not  can  be  determined  intuitively  from 
the  observed  pattern  of  concentration  distribution.  One  should  also  be 
concerned  with  the  effect  of  the  bottom  boundary  so  that  the  diffusion 
in  question  may  be  classified  as  two-dimensional  or  three- 
dimensional.  For  this  purpose,  a  sharp  thermocline  is  regarded  as  a 
boundary.  The  presence  or  absence  of  constraints  in  the  vertical 
direction  does  not  create  a  drastic  'vange  in  horizontal  diffusion  but 
have  a  significant  influence  on  the  time  behavior  of  the  maximum 
concentration  of  patch. 

3)  The  horizontal  distribution  of  dye  concentration  should  be 
observed  in  such  a  manner  so  that  the  horizontal  variance  associated 
with  the  dye  distribution  can  be  computed  directly.  Thus  aerial 
photographs  of  the  dye  patch  are  not  considered  useful  for  this 
purpose,  although  some  progress  has  been  made  in  determining  the 
concentration  distribution  from  dye  photos  (Ichiye  and  Plutchak,  1966). 
In  observing  the  time  change  of  patch  size  in  photographs,  Gifford 
(1957)  proposed  an  indirect  method  to  estimate  the  variance  in  the 
case  of  atmospheric  diffusion.  Ichiye  (1959)  and  Okubo  (1962) 
applied  this  technique,  seemingly  with  some  success,  to  horizontal 
diffusion  of  dye  in  the  sea.  However,  Gifford's  technique  is  not 
employed  in  this  report,  not  only  because  the  direct  computation  of 
the  variance  is  more  advantageous  but  also  because  the  variance 
estimated  by  means  of  the  Gifford  technique  depends,  to  some  extent, 
on  the  choice  of  the  assumed  distribution  for  dye  concentration. 

4)  The  concentration  distribution  should  agree  with  a  "mass  or 
material  balance",  by  which  is  meant  that  the  amount  of  dye 
remaining  within  a  patch  must  be  comparable  to  the  total  amount  of 
dye  released,  when  .he  loss  of  dye  due  to  photochemical  decay, 
physical  adsorption,  etc.  ,  if  appreciable,  is  corrected.  Pritchard 
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and  Carpenter  (I960)  in  their  study  of  diffusion  in  Chesapeake  Bay, 
reported  that  mass  balances  of  80  to  90%  were  realized  using 
Rhodamine  B  as  a  tracer.  This  case  is,  however,  a  little  ideal. 

Usually  mass  balances  as  low  as  80%  or  as  high  as  150%  may  be 
taken  as  permissible  especially  in  open  ocean  experiments  and  in  the 
later  stage  of  diffusion  characterized  by  very  low  concentrations  of 
dye  within  a  patch.  Quite  often  the  mass  balance  is  neither  reported  ‘j 

nor  computable  simply  because  of  lack  of  information  on  the  vertical 
distribution  of  dye.  In  this  case  an  appropriate  value  of  the  mixed 
layer  within  which  dye  is  assumed  to  be  distributed  uniformly  is  taken 
on  the  basis  of  some  indirect  information,  say  the  stratification  of  the 
water  column.  The  amount  of  dye  in  a  patch  is  then  computed  from 
the  horizontal  distribution  of  dye  multiplied  by  the  depth  of  the  mixed 

layer.  Another  way  of  examining  the  mass  balance  represents  an  i 

i 

indirect  approach:  take  the  observed  horizontal  distribution  of  dye 
and  compute  the  depth  of  mixed  layer  within  which  100%  of  material 

r  1 

is  supposed  to  remain.  If  this  fictitious  depth  turns  out  to  be 
reasonable  from  the  standpoint  of  the  knowledge  of  vertical  diffusion 
in  the  sea,  the  stratification  of  the  water  column,  the  total  depth  of 
water,  etc.  ,  we  may  regard  the  experiment  as  "good"  with  respect  to 
a  mass  balance. 

i 

i 

3.  Computation  of  diffusion  characteristics.  A  diffusion  diagram  I 

consists  of  a  plot  of  a  characteristic  property  of  diffusion  against  an 

! 

another  such  property.  A  set  of  the  most  natural  choices  for  1 

characteristic  properties  is  the  variance  of  the  horizontal  distribution*  l 

.  i 

of  material  versus  the  diffusion  time,  i.  e.  ,  the  time  elapsed  since  the  j 


* 

The  maximum  concentration,  for  example,  may  not  be  a  good 
measure  of  diffusion,  not  only  because  the  observation  of  the 
maximum  concentration  involves  a  great  deal  of  uncertainty  but  also 
because  the  peak  concentration  is  very  sensitive  to  the  presence  or 
absence  of  a  horizontal  boundary. 
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introduction  of  material.  Usually  the  origin  of  the  time  (the  "  release 
time")  is  defined  as  the  midpoint  of  the  time  of  release,  and  the 
diffusion  time  is  the  difference  between  the  midpoint  of  the  time  spent 
in  observing  the  concentration  distribution  and  the  release  time.  The 
part  of  the  concentration  distribution  observed  prior  to  a  diffusion 
time  would  be  diffused,  to  some  extent,  if  one  observed  it 
instantaneously  at  the  diffusion  time.  On  the  contrary,  the  other  part 
of  the  distribution  observed  after  the  diffusion  time  would  be  less 
diffused  if  observed  instantaneously  at  the  diffusion  time.  Thus  the 
net  effect  is  compensatory,  so  that  a  more  or  less  instantaneous  value 
of  the  variance  may  be  obtained  with  respect  to  the  diffusion  time. 

The  variance  of  the  horizontal  distribution  is  certainly  a  suitable 
measure  of  the  spread  of  the  contaminant.  Horizontal  diffusion  in  the 
sea  can  not  be  described  adequately  by  a  Fickian-type  diffusion  with  a 
constant  coefficient  of  diffusion.  Thus  the  theory  of  oceanic  diffusion 
has  been  directed  toward  non-Ficklan  diffusion;  as  a  matter  of  fact, 
the  apparent  diffualvity  increases  with  the  diffusion  time  or  with  the 
size  of  the  diffusive  substance  (OkuDO,  1962). 

In  order  to  estimate  the  variance  directly,  one  must  obtain  a 
fairly  good  pattern  of  the  horizontal  concentration  distribution  at  a 
particular  depth.  An  example  of  such  a  concentration  distribution  is 
shown  in  Fig.  1,  where  we  also  indicate  the  ship  tracks  along  which 
continuous  recordings  of  concentration  are  made.  Since  the 
horizontal  distribution  of  substance  is  usually  asymmetric  in  the 
sense  that  a  characteristic  length  is  larger  in  one  direction  or  another, 
two  variances,  along  the  long  and  short  axes,  say,  are  necessary  to 
properly  describe  the  rate  of  dispersion.  However,  conversion  to  a 
radially  symmetric  equivalent  distribution  has  been  widely  used  in 
reporting  the  results  of  dye  experiments. 

The  concept  involved  in  a  radially  symmetric  patch  is  as 
follows.  At  a  diffusion  time  t,  the  shape  of  the  isolines  of  any 
concentration  will  be  very  irregular  in  any  individual  case;  the  patch  of 
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I  substance  is  often  elongated.  We  shall  now  carry  out  (conceptually!) 

|  an  infinite  number  of  releases  under  identical  oceanographic 

i 

conditions,  using  in  each  case  the  same  amount  of  substance  released. 
Thus,  if  we  take  an  infinite  number  of  distributions  each  observed  at 
.  the  same  diffusion  time,  if  we  superpose  them  in  such  a  way  that  the 

center  of  mass  coincides  with  each  other,  and  if  we  average  over  all 
the  superposed  distributions,  we  would  then  expect  a  radially 
symmetric  distribution  of  substance  about  the  center  of  mass, 
provided  that  the  oceanic  conditions,  e.  g.  ,  the  shears  in  the  mean 
current,  are  isotropic. 

Normally  identical  oceanic  conditions  do  not  necessarily  imply 
rotational  symmetry  in  the  shears.  A  case  in  which  the  shears  in 
the  mean  flow  are  likely  to  occur  predominantly  in  one  direction  is 
seen  in  the  nearshore  region.  There  almost  always  a  patch  of  dye 
tends  to  elongate  in  the  direction  parallel  to  the  coast.  Thus  the 
ensemble  average  of  distribution  would  be  more  or  less  elongated 
rather  than  radially  symmetrical.  A  radially  symmetrical 
distribution  is,  however,  the  simplest  mathematical  model  to  deal 
with.  Thus  the  variance  associated  with  a  radially  symmetric 
distribution  can  be  related  through  a  measure  of  the  eccentricity  to 
the  variances  for  a  two-dimensional  elliptical  distribution  which 
itself  is  a  proper  mathematical  model  for  an  elongated  patch  of  dye 
(Okubo,  1965).  The  Appendix  summarises  the  relationships  between 
various  variances  associated  with  particular  mathematical  models  for 
the  two-dimensional  distribution  of  substance. 

Intuition  tells  us  that  the  area  enclosed  by  a  line  of  constant 
concentration  would  be  nearly  the  same  for  all  identical  releases  in 
spite  of  the  very  irregular  distribution  of  concentration  in  individual 
releases.  As  a  measure  of  the  area  enclosed  by  a  particular 
concentration,  we  define  th«  radius  r  of  a  circle  of  equal  area.  A 
radially  symmetrical  distribution,  S  (t,  r^),  thus  characterized  by 
the  equivalent  radius  and  the  diffusion  time  would  represent 
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approximately  the  radially  symmetric  distribution  in  the  ensemble 
average  mentioned  previously. 

Having  obtained  a  sufficient  number  of  data  on  r  vs.  S  for  a 

C 

2 

patch  of  dye,  one  may  compute  the  variance  afo  by  the  definition. 

00  00 

°ro2(t)  E  J  re2  8(t're)  2nTe  ‘V  J  S  2nre  ^e  ’  (2) 

o  o 


Frequently  the  spatial  distribution  of  S  is  well  described  by  a  two* 

dimensional  Gaussian  distribution,  i.  e.  , 

2 

r 

e 


S(t,re) 


-M/P—  e 

M  Z(t) 

TO 


) 


(3) 


where  M/D  denotes  the  total  mass  of  dye  per  unit  depth  of  mixed 

layer.  In  this  case  ar<j2  can  be  estimated  from  a  plot  of  log  S  vs. 

r^2.  Fig.  2  shows  such  an  example. 

In  another  case  S(t,  r  )  can  be  described  by  a  generalized 

G 

exponential  distribution: 


S(t,re) 


r  m 

S(t,o)  e  r° 


(4) 


whore  S(t,  o)  is  the  peak  concentration  and  m  is  a  positive  number. 


Thus  a  plot  of  log  9.1  vs.  r  will  determine  the  value  of 

2  Svt,rgl  e 

oro  as  well  as  m.  An  example  is  seen  in  Fig.  3.  In  any  case,  we 

2 

can  calculate  <j  directly  from  definition  (2). 
re 

2 

The  accuracy  of  the  estimate  of  aro  depends  primarily  on 
uncertainties  in  the  equivalent  radii.  "Good"  data  by  which  is  meant 

that  a  fairly  precise  distribution  of  concentration  is  observed,  should 

2 

provide  an  estimate  of  orc  within  ±  5%  of  accuracy.  Completely 


Fig.  3.  Fit  of  diffusion  data  to  a  generalized  exponential 
distribution. 
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independent  evaluations  of  the  equivalent  radii  by  the  use  of  the  same 

2 

raw  data  reveals  that  the  relative  difference  in  r^  is  a  few  percents. 

The  early  stages  of  dispersion  are  difficult  to  sample  since 
the  patch  has  small  areal  extent  and  any  significant  disturbance 
of  tht  patch  by  a  sampling  vessel  must  be  avoided.  An  example  shows 
clearly  how  a  small  disturbance  of  unknown  cause  created  a  trigger 
action  on  the  diffusion  of  patch.  Figs.  4a  through  4d  are  serial 
photographs  of  a  dye  patch  at  early  time.  It  may  be  seen  from  the 
figures  that  at  time  between  4m  51*  and  6m  47*  after  release,  the 
leading  edge  of  the  patch  turned  180*  while  lagging  behind  and  moving 
relative  to  the  trianglar  drogue  with  a  length  of  7.  3  meters.  The 
existence  of  small-scale  eddying  motion  is  evident.  By  9m  08*  after 
release  the  leading  edge  had  rejoined  the  main  body  of  the  patch 
resulting  in  an  elliptical  shape.  After  lliri  408  after  release,  the 
mixing  of  dye  had  taken  place  internally  so  that  a  fairly  large  patch  of 
50  meters  in  scale  appeared  (for  details  see  Carter,  1967).  A  simple 
calculation  baaed  on  Joseph  and  Sendner  theory  with  P  =  1  cm/s  shows 
that  it  would  have  taken  more  than  an  hour  for  the  initial  patch  to  grow 
into  a  patch  of  50  meters  in  diameter  without  the  effect  of  the 
disturbance. 

With  a  view  toward  minimizing  the  disturbances  due  to  a 
sampling  vessel  yet  still  determining  the  variances  of  the 
concentration  distribution,  Foxworthy  et  al  (1966)  carried  out  dye 
release  experiments  in  which  a  small  sampling  vessel  made  single 
crossings  of  dye  patches  of  small  size  down  to  30  meters.  The 
crossings  were  aimed  so  as  to  lie  along  the  direction  of  elongation  of 
a  patch  and  the  perpendicular  direction  so  that  the  variances  along 
the  principal  axes  of  the  concentration  distribution  could  be 
calculated.  The  observed  concentration  distributions  have  been 
fitted  with  a  Gaussian  curve  to  compute  the  variances.  The  fit  was 
found  to  be  reasonably  good  especially  during  the  early  phase  of 
diffusion. 


.  6m  47s  after  release 


9m  03s  after  release 


ft 


11m  40s  after  release 

Fig.  4  Aerial  photographs 


study  in  Manokin  River 
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Since  the  two  crossings  along  the  principal  axes  were  not 
carried  out  at  the  same  time,  an  empirical  relationship  between  the 
estimated  variance  against  diffusion  time  was  determined  for  each 
a  particular  axis.  A  log-log  plot  of  the  variance  against  time 
determines  a  straight  line  which  provides  the  empirical  law  expressed 


2  ,  .  n 

°x  “kt 


say  for  the  variance  along  the  major  principal  axis.  It  can  be  seen 
(Appendix)  that,  if  crx2  and  oy2  denote  the  one -dimensional 

variances  along  the  principal  axes  of  a  two-dimensional  normal 

2 

distribution,  the  variance  orc  for  a  radially  symmetrical 
distribution  obtained  by  taking  the  equivalent  radius  is  given  by 


2  a  a 
x  y 


Actually  dye  patches  elongate  in  one  direction  or  another,  and  the 

distribution  of  concentration  may  be  represented  reasonably  well  by 

a  two-dimensional  Gaussian  distribution.  Thus,  from  the  empirical 

laws  of  o 2  and  o  2  with  time,  we  can  compute  o  2  for 
x  y  rc 

specified  times  of  diffusion. 

All  other  dye  experiments  herein  quoted  consisted  of 
measuring  the  horizontal  distribution  by  a  proper  network  of 
crossings  in  a  dye  patch.  When  the  actual  pattern  of  dye  distribution 
is  shown  in  reports  (e.  g.  ,  Joseph,  Sendner  and  Weidemann,  1964), 

one  may  compute  the  mean  variance  a  2  ,  the  individual  variances 

2  2 

a_  and  o„  ,  the  degree  of  elongation  a  /a.  etc.  ,  in 


addition  to  the  variance  for  a  radially  symmetrical  distribution. 

,  <J 

An  apparent  diffusivity  K&  defined  by  i  rc  and  a 


scale  of  diffusion  defined  by  3  o 


are  also  computed;  the  scale  of 
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diffusion  is  arbitrarily  taken  as  the  diameter  of  a  radially  symmetric 
patch  within  which  95%  of  substance  remains  (see  Appendix).  The 
results  of  the  computations  are  summarized  in  Table  1. 

4,  Diffusion  diagrams.  By  plotting  the  value  of  the  variance  a  2 

.  . .  "  ■  a  r  c 

against  the  time  of  diffusion,  we  obtain  a  basic  diffusion- diagram 

(Fig.  5).  The  variance  ranges  from  10^  cm^  up  to  10^  cm*",  while 

the  diffusion  time  ranges  from  2  hours  to  nearly  a  month.  For 

reference  purposes,  we  include  the  diffusion  data  in  the  Manokin 

River  (a  narrow  tidal  tributary  to  Chesapeake  Bay)  and  the  Banana 

River  (s  shallow  basin).  Horizontal  diffusion  in  these  regimes 

should  be  regarded  as  two-dimensional.  Though  the  points  in  the 

diagram  (Fig.  5)  scatter  somewhat,  an  obvious  trend  is  recognized 

that  the  variance  grows  much  faster  in  the  power  of  t  than  a  linear 

increase  with  time;  a  Fickian  diffusion  characterized  by  a  constant 

diffusivity  would  result  in  the  variance  which  grows  linearly  with 

diffusion  time.  Thus  horizontal  diffusion  in  the  sea  is  a  process  in 

which  the  apparent  diffusivity  would  increase  with  the  time  of 

diffusion  or  with  the  scale  of  diffusion.  A  linear  fit  to  all  the  data 

points  in  Fig.  5  gives  the  relationship  that 

a  2  =  0.0106  •  t~‘54  (a  Z\  cm2,  t:  sec)  .  (5) 

rc  rc  ’ 


Using  diffusion  data  prior  to  1961,  Okubo  (1962)  obtained 

2 

empirically  the  similar  relationship  between  a ^  and  t: 

o  2  =  0.006  •  t?,'>  (a  2 :  cm2,  t;  sec) 

rc  rc 


Although  the  quality  of  the  old  diffusion  data  is  poor  in 
comparison  with  that  of  the  new  date  used  in  this  report,  the  general 


Table  1  Computation  of  diffusion  characteristics 
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Degree  of  elongation 


100  km 


Fig,  5.  Variance,  a 

r  o  ' 


v».  diffueion  time,  (new  data) 
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conclusion  remains  unchanged.  Although  the  new  data  may  be  fitted 
locally  in  time  to  the  power  3.  0,  the  overall  slope  of  the  log  o^-log  t 
line  is  noticeably  smaller  than  3.  0;  the  third  power  of  t  for  the 
variance  has  been  predicted  by  the  similarity  theory  of  turbulence 
(see  Section  5). 

In  Fig.  6  we  plot  all  the  data,  old  and  new  together,  in  one 
diagram.  It  should  be  noticed  that  the  observed  variance  at  a 
diffusion  time  varies  as  much  as  an  order  of  magnitude.  The 
variability  in  the  variance  seems  to  be  more  noticeable  at  smaller 
diffusion  times.  For  large  times  of  diffusion,  a  patch  of  diffusing 
substance  itself  takes  average  conditions;  a  variety  of  wind  and 
current  situations  is  encountered  and  the  average  effect  does  not 
differ  much  from  one  experiment  to  another. 

Fig.  7  shows  the  relationship  between  the  apparent  diffuslvity 
and  the  scale  of  diffusion.  A  linear  fit  to  the  new  data  points  gives 
the  relationship  that 


Ka  =  0.0103  *  i1*15  (Kft:  cm2/*,  /:  cm)  .  (6) 


The  exponent  of  L  is  definitely  smaller  than  4/3;  the  4/3  law  of 
oceanic  diffusion  has  been  proposed  by  several  investigators  (see 
Bowden,  1962),  and  explained  by  them  with  the  similarity  theory  of 
turbulence.  So  far  as  dye  diffusion  experiments  are  concerned,  how¬ 
ever,  the  apparent  diffusivity  does  not  obey  the  4/3  law,  although 
there  remains  the  possibility  that  the  4/3  law  may  be  applied  locally 
in  scale  to  the  data.  Okubo  (1962)  came  to  the  same  conclusion  by 
the  use  of  the  old  data;  the  old  relationship  was  given  by 
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Ko  ( cm2  /g#c) 
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Fig.  7.  Apparent  diffusivity,  K&,  vs.  scale  of  diffusion,  l  =  3  apc 
(new  data). 
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.  d  e 

where  K  was  defined  by  •=  — — 
r  2  dt 


and  t  by  4  a 


In  Fig.  8  we  plot  all  the  data,  old  and  new,  together  in  one  diagram 
of  the  K  -  l  relationship.  The  diagrams  (Figs.  5  and  7)  can  be 

2L 

used  for  ari  empirical  approach  to  predict  the  horizontal  spread  of 
substance  from  an  instantaneous,  small  source,  within  a  certain 
range  of  variability  in  the  <  '  -iiacteristics.  An  example  of  predictions 
by  the  use  of  the  relations  (5)  and  (6)  is  given  in  Table  2. 


5.  Discussion.  The  new  diffusion  diagrams  herein  constructed  using 
the  data  of  recent  dye  experiments  provide  not  only  a  practical  means 
of  predicting  the  characteristics  of  oceanic  diffusion  but  also  a 
certain  clue  concerning  the  spectrum  of  oceanic  turbulence.  As  for 
the  prediction  of  horizontal  diffusion  from  an  instantaneous  point 
source,  Figs.  5  and  7  can  be  used  as  a  general  guide.  Irrespective 
of  the  detailed  oceanographic  conditions,  the  horizortal  diffusion 
exhibits  a  general  trend;  the  variance  increases  widi  time  at  a  power 
be*  sen  2  and  3.  1.  is,  however,  to  be  remembered  that  the 

variance  may  vary.probably  depending  on  the  oceanographical 
conditions,  by  an  order  of  magnitude  or  more  for  the  same  diffusion 
time.  Only  a  limited  number  of  dye  studies  have  been  made  in  which 
the  horizontal  spread  of  substance  is  compared  with  the  environmental 
factors,  such  as  the  stability  of  water  column  (Foxworthy  et  al,  1966, 
Kullenberg,  1967). 

Foxworthy  et  al  (1966)  examined  the  effect  of  water  column 
stability  and  wind  speed  on  the  variances  in  the  directions  of  the 
principal  axes  of  diffusing  patch.  At  any  given  average  value  of  wind 
speed  ranging  from  2  to  14  kt,  lower  values  of  dilution  (a  measure  of 
the  rate  of  dispersion)  are  associated  with  higher  values  of  average 
stability.  This  effect  is  controlled  primarily  by  the  suppression  of 
the  vertical  diffusion  by  increasing  stability.  The  effect  of  stability 
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Table  2  An  example  of  predictions  of  diffusion  characteristics 


a)  Variance 

and  standard 

deviation  of  a  radially  symmetric  patch 

vs  diffusion  time 

t(sec) 

2  \ 

<Jpc  (cm  ) 

a  (cm) 
rc 

3.6  xlO3 

(l  hour) 

2.3  xlO6 

1.5xl03 

io4 

2-5  xlO7 

5.0xl03 

3.6  xlO4 

(10  hoars) 

4.9  xlO8 

2.2xl04 

e.64xio4 

(1  day) 

3-6  xlO9 

6.0xl04 

105 

5-4  xlO9 

7.4xl04 

6.03xJ.o5 

(1  week) 

3-6  xlO11 

6.Qxl05 

IO6 

1.2  xlO12 

l.lxlO6 

2. 59x10 s 

(l  month) 

1.2  XlO13 

3-3xl06 

b)  Apparent  diffusivity 

vs  scale  of  diffusion 

i(cm) 

Ka(cmVsec) 

io4 

4.LxlO? 

5  x  IO4 

2.6x103 

105 

5.8xl03 

5  x  io5 

3.7xl04 

io6 

8.2X104 

5  x  IO6 

5-3xl05 

IO7 

1 . 2xl06 
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on  horizontal  dispersion  is  inconclusive,  however.  Thus,  in  the  range 
of  2  to  8  kt  wind,  an  increase  in  o^  with  decreasing  stability  is 
suggested,  but  there  are  insufficient  data  on  which  to  base  a  firm 
conclusion.  Similar  stability  appears  to  have  little  or  no  definite 
effect  on  o^  .  Although  the  data  scatter  considerably  on  a  plot  of 
the  variances  and  wind  speed,  there  is  indicated  a  trend  of  increasing 
a  2  with  increasing  wind  speed,  whereas  there  appears  to  be  no 
definite  correlation  between  o  and  wind  speed.  The  wind  speed 
of  our  concern  may  be  regarded  as  a  measure  of  a  characteristic 
current  speed  in  the  water  column  within  which  diffusion  takes  place. 

It  has  recently  become  apparent  that  velocity  shear  can  play  an 
important  role  in  the  horizontal  spread  of  contaminants  in  the  sea. 

The  basic  result  of  the  effect  of  current  shear  is  that  an  effective 
longitudinal  dispersion  is  produced  by  the  combination  of  the  gradient 
of  mean  velocity  and  turbulent  mixing  in  the  same  direction.  In  the 

case  of  a  uniform  vertical  shear  in  the  horizontal  current,  the 

2  2 

longitudinal  and  lateral  variances,  ax  and  cy  ,  can  be 
expressed  as  (Carter  and  Okubo,  1965): 


a  = 


i- 


2  A  t’ 
z  z 


2  A  t 
y 


(7) 


where  n  :  the  constant  vertical  shear,  A  :  the  vertical  eddy 
z  z 

diffusivity.  Ay:  the  lateral  eddy  diffusivlty. 

It  may  be  speculated  that  the  lateral  eddy  diffusivity  would  be 
relatively  unaffected  by  stability.  That  is,  would  be 

relatively  uncorrelated  to  stability.  On  the  other  hand,  the  tendency 
of  increasing  a  with  increasing  wind  speed,  i.  e.  ,  increasing 
shear,  is  evident  from  the  expression  (7),  and  the  decrease  of  cr^2 
with  increasing  stability  can  also  be  seen  because  of  the  suppression 
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of  the  vertical  eddy  diffuaivity  by  increasing  stability.  Thus  the 
simple  shear-diffusion  model  can  interpret,  at  least  qualitatively, 
the  experimental  results  discovered  by  Foxworthy  et  al.  The  shear 
effect  on  horizontal  diffusion  implies  that  the  variance  for  a  radially 
symmetrical  distribution  is  closely  related  to  the  current  shear  and 
stability  aside  from  the  time  dependence.  The  variability  of  the 
variance  at  a  fixed  diffusion  time  is  thus  nothing  but  a  reflection  of 
the  change  in  fhe  environmental  conditions.  A  future  study  should 
include  these  environmental  factors  as  parameters  in  the  discus3;or. 
of  the  diffusion  diagram. 

As  mentioned  previously,  a  radially  symmetrical  distribution 
may  not  represent  properly  the  actual  pattern  of  oceanic  diffusion. 

In  addition,  the  value  of  the  horizontal  variance  is  always  smaller 
than  it  should  be  when  we  convert  the  actual  distribution  of 
concentration  S(x,y,t),  into  a  radially  symmetrical  one  characterised 
by  an  equivalent  radius  (see  Appendix).  In  other  words,  we 
underestimate  the  horizontal  variance  by  an  arbitrary  amount  when 
taking  a  radially  symmetrical  distribution.  The  correction  factor 
/a  ^  depends  on  the  degree  of  asymmetry  in  the  actual 

distributior  (Appendix).  The  more  elongated  the  actual  patch  of 
substance,  the  larger  the  correction  factor  should  be. 

We  have  observed  quite  often  that  the  patch  tends  to  elongate 
more  and  irm-e  s  time  goes  on  (see  a  paper  by  Joseph,  Sendner  and 
Weidemann,  1964).  Thus  the  correction  factor  would  grow  as  the 
diffusion  time  grows.  This  raises  an  important  question;  "How 

2  2  I’ 

would  a  behave  with  time  if  we  plot  <3  ,  instead  of  0  , 

r  r  re 

against  t0"  At  present  not  much  can  be  answered  because  some  of 

the  large  scale  diffusion  data  provide  information  only  on  the 

2 

equivalent  radius.  Fig.  9  shows  a  plot  of  a  vs.  t.  It  suggests 

3  2 

more  clearly  than  Fig.  5  that  the  t  law  of  or  should  be  valid 
locally  in  time. 

The  turbulent  diffusion  of  a  patch  of  dye  is  a  problem  in  relative 
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diffusion.  A  theoretical  treatment  of  this  problem  has  been  given  by 
Batchelor  (1952).  With  the  aid  of  dimensional  arguments.  Batchelor 
finds  a  regime  of  relative  diffusion  if  the  initial  size  of  the  patch  is 
infinitesimal.  This  regime  is  characterized  by  the  following 


CEt 


(8) 


where  C  is  a  constant  of  order  unity  and  Els  the  rate  of  energy 
dissipation  per  unit  mass  through  turbulence.  It  is  easy  to  derive 
from  (8)  that  the  apparent  diffusivity  grows  as  the  4/3  power  of 
patch  size: 


K 


El/3  o4/3 
r 


(9) 


This  is  Richardson's  (1926)  well-known  law  of  relative  diffusion,  the 
physical  meaning  of  which  is  that  relative  diffusion  is  an  accelerating 
process,  the  rate  of  growth  increasing  with  the  size  of  the  patch. 

The  basic  concept  involved  in  the  relationship  (8)  and  also  in  the 
Richardson's  law  is  that  the  eddies  responsible  for  the  horizontal 
spread  of  substance  lie  in  the  "inertial  subrange"  (Kolmogorov,  1941). 
These  eddies  receive  their  energy  from  the  larger  eddies  and  pass  it 
to  eddies  smaller  than  themselves.  No  direct  energy  is  supplied  to 
those  eddies  from  external  sources  and  the  energy  dissipation  due  to 
viscosity  is  not  significant.  Their  properties  thus  depend  only  on  the 
rate  of  energy  transfer  which  must  be  equal  to  the  rate  of  energy 
dissipation  if  the  energy  of  the  eddies  remains  stationary. 

Under  natural  conditions,  however,  there  is  often  the 
possibility  of  energy  being  added  directly  to  such  eddies  by  wind 
systems,  tidal  currents,  or,  on  a  still  smaller  scale,  by  waves 
(Stommel,  1949,  Ozmidov,  1965).  These  excited  eddies  could  not 
then  be  strictly  in  the  inertial  subrange,  so  that  results  derived  from 
that  concept  would  not  necessarily  be  applicable.  Ozmidov  (1965), 
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nevertheless,  suggests  the  existence  of  locally  Inertial  subranges 
separated  by  two  scales  of  eddies  where  an  influx  of  external  energy 
takes  place.  Fig.  10  shows  schematically  the  spectrum  conjectured 
by  Ozmidov.  This,  in  turn,  suggests  that  the  law  (8)  or  the 
Richardson's  law  may  be  used  locally  in  time  or  scale  with  a  different 
value  of  the  rate  of  turbulence  energy  transfer.  As  the  scale  of 
diffusion  Increases,  the  smaller  the  value  of  this  parameter  one 
should  expect,  since  the  local  supplies  of  energy  tend  to  be 
transferred,  through  the  nonlinear  Interactions,  from  the  larger 
eddies  to  the  smaller  eddies  (for  some  exceptional  cases,  see  Webster 
1961).  As  a  consequence,  the  overall  rate  of  growth  of  the  variance 
with  time  should  be  slower  than  the  cube  power  law  and  the  overall 
growth  rate  of  the  apparent  diffusivity  with  the  scale  should  be 
smaller  than  the  4/3  law. 

On  the  basis  of  this  concept,  we  fit  the  law  (8)  locally  to  the  data 

points  as  shown  in  Fig.  6,  from  which  the  values  of  the  rate  of  energy 

transfer  in  the  local  regimes  are  computed,  assuming  the  numerical 

constant  C  a  1.  Similar  fit  to  the  data  of  K  vs.  £  is  shown  in  Fig.  8. 

a 

The  result  of  computation  of  E  is  summarized  by 
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A  clear  jump  in  the  value  of  E  occurs  at  a  time  of  a  half  day  to 
a  day.  This  suggests  a  local  energy  supply  by  tidal  motions, 
amounting  to  9.  4  x  10  ^  cm^/sec^  per  unit  mass  of  sea  water. 
According  to  Jeffrey  (1952)  the  average  rate  of  dissipation  of  tidal 
energy  is  estimated  at  1.4  x  ergs/sec.  If  this  were  distributed 

uniformly  throughout  the  ocean  body,  the  average  rate  of  dissipation 


of  tidal  energy  per  unit  mass  of  water  would  be  1  x  10'^  cm^/sec^, 

.5 

which  is  an  order  of  magnitude  smaller  than  the  value  of  9.4  x  10 
cm^/sec^.  However,  it  is  believed  (Knauss,  1956)  that  the  most  of 
the  tidal  energy  is  dissipated  in  shallow  seas;  thus  the  local  energy 
supply  by  tidal  motions  would  be  much  higher  than  10  ^  cm^/sec^  in 
shallow  waters. 

There  appears  to  exist  another  jump  in  E  at  a  time  scale  of  a 

few  to  several  days.  This  is,  however,  inconclusive  not  only  because 

of  the  fact  that  the  data  points  for  the  large-Bcale  diffusion  have  not 

been  corrected  to  the  mean  variance,  a  ,  but  also  because  the 

r 

quality  of  those  points  which  cover  a  scale  of  more  than  100  km  is 
found  to  be  very  poor. 

6.  Final  Remarks. 

1)  The  diffusion  diagrams  herein  presented  are  made  exc.  jively 
from  the  data  of  diffusion  experiments  in  the  surface  layer  of  the  sea. 
Very  little  has  been  known  about  the  horizontal  diffusion  in  deep  water 
or  in  the  thermocline  simply  because  of  the  technical  difficulties  in 
carrying  out  dye  release  experiments  there.  Doubtless  those 
experiments  would  provide  important  clues  to  our  understanding  of 
large-scale  mixing  in  the  sea. 

2)  In  order  to  provide  adequate  information  on  the  study  of 
oceanic  diffusion,  a  particular  manner  of  data  reporting  must  be  used. 
A  good  example  of  data  presentation  is  seen  in  papers  by  Joseph, 
Sendner  and  Weidemann  (1964)  and  by  Carter  and  Okubo  (1965).  The 
following  information  should  be  considered  as  "necessary",  whereas 
some  of  them  have  often  been  missing  in  reports; 

i)  Vertical  distribution  of  concentration  as  a  measure 
of  depth  of  mixing, 

ii)  Velocity  field,  especially  current  shears  in  the  scale  of 
diffusion;  a  minimum  requirement  is  the  mean  velocity  of  the 
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center  of  mass  of  patch  during  the  interval  of  successive 
observations , 

iii)  Initial  conditions  of  release:  the  duration  of  release,  the 
depth  of  release,  the  initial  size  of  patch,  etc.  , 

iv)  Shape  information:  the  degree  of  elongation,  the  direc¬ 
tion  of  elongation  with  respect  to  the  mean  flow,  etc.  ;  the 
pattern  of  the  horizontal  distribution  of  concentration  should 
be  shown  whenever  available, 

v)  Background  oceanographic  conditions:  the  stability  of 
water,  etc.  , 

vi)  Wind  data,  sea  state,  etc. 
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APPENDIX 

Characteristics  of  dispersion  associated  with  a  two-dimensional 
distribution 


Consider  a  two-dimensional  distribution  of  concentration  from 
an  instantaneous  release  of  a  unit  amount, 


S  =  S(t,x,y):  ff  SAxdy  =  1  for  all  t  2  0  .  (A-l) 

-00 


Thus  the  mean  value  of  any  property  f(x,  y,  t)  associated  with  S  is 
defined  by 


f 


i  r* 

fSdxdy  . 


-00 


(A -Z) 


The  origin  of  the  orthogonal  coordinate  system  (x,  y)  with  an 
arbitrary  orientation  is  taken  at  the  center  of  mass  of  the  distribution 
S,  so  that  x  =  y  =  0.  The  mean  square  distance  from  the  center 
of  mass,  i.  e.  ,  the  variance,  defining  by 


o  ?  s  r2  -  f f  (x2  +  y2)  Sdxdy  , 
r  dc 


(A-3) 


thus  can  be  expressed  as  the  sum  of  the  mean  square  distances  from 
the  center  of  mass  along  the  two  orthogonal  coordinates; 


2  2  .  2 
a  =0  +0 

r  x  y 
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Let  us  compute  the  mean  square  separations  between  a  pair  of 
subetance  particles  in  the  patch,  l 2  : 


ta(*v 
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|2  =  JJW  kJ2  ^  » 


(A-5) 


where  and  are  the  position  vectors  from  the  origin  to 

two  particles,  respectively,  £  “  ,  and 

is  the  joint  probability  density  function  for  the  two  particles  (see  Fig. 
A).  The  integral  (A-5)  extends  over  the  whole  space  of  ^  and  r.. 

Since  we  can  take  two  particles  each  independently  of  the  other, 
the  joint  probability  is  expressed  as  the  product  of  the  density 
functions  for  each  particle,  that  is 


toiSjjJCg)  =  3(t>c1)  8(1^) 


(A-6) 


Substitution  (A-6)  into  (A-5)  and  calculating,  we  obtain 


f  =  J  X22  S(t,r2)dr2  J  3(6,^)*^ 


2  PP  r  •  r  S(t,r  )  S(t,r  )  dr  dr. 
JJ  »v2  M  'V’l 


+  J  »(«<*,>  %  J  2t, 


(A-v) 


£  2  +  j^2  (since  =  0)  =  r2  +  r2  =  2  o/  . 
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We  thua  conclude  that  the  mean  square  separations  between  a  pair  of 
particles  in  the  patch  is  equal  to  the  twice  of  the  mean  square 
distances  of  a  particle  from  the  center  of  mass  of  the  patch. 

We  shall  now  specify  the  distribution  by  a  two-dimensional 
Gaussian  distribution: 


S(t,x,y)  = 


l-X  2o. 


2«o1(t)o2(t)(l  -  X^t))" 


(A-8) 


where  a^Ct)  and  are  the  variances  in  the  directions  of  x 

and  y,  respectively,  and  X  is  a  coefficient  of  correlation: 


(A-9) 


Equation  (A-8)  describes  that  the  concentration  contours  are  a  set  of 
ellipses  with  respect  to  the  center  of  mass. 

Rotating  the  coordinate  system  (x,  y)  by  a  certain  angle,  we  can 
express  S  in  terms  of  the  principal  axes  (X,  Y). 


S(t,X,Y) 
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(A- 10) 
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where  and  a  are  the  variances  in  the  major  and  minor 

principal  axes,  respectively.  Obviously, 
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We  now  convert  the  asymmetric  distribution  (A- 10)  into  a 

radially  symmetric  one  by  taking  the  equivalent  radius,  r  .  Since 

© 

the  area  enclosed  by  a  concentration  must  be  equal  for  both 
distributions,  we 
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Hence,  (A- 10)  is  reduced  to 
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Define  the  variance  for  the  radially  symmetric  distribution  by 
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We  then  obtain 
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Therefore, 
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In  other  words,  the  variance  for  the  radially  symmetric  distribution 
is  equal  to  the  twice  of  the  geometric  average  of  the  one -dimensional 
variances  in  the  principal  axes  for  the  original  two-dimeneional 
Gaussian  distribution. 

It  should  be  noted  that,  in  general, 

°rc?^2Vy  ’  (A- 15) 

where  x  and  y  are  orthogonal  coordinates  in  an  arbitrary 
orientation  with  respect  to  an  elliptical  patch  of  uubstance. 
i  It  is  easily  seen  that  the  variance  for  an  elliptical  distribution 

never  amounts  to  less  than  that  of  the  corresponding  radially 
symmetrical  distribution,  that  is 

2  2  2  2  /  \  2  ^  ^ 
a  -a  =  ov  +  o„  -  2a  a  =  (o  -  a)  SO  . 
r  ro  X  Y  X  Y  X  Y 

Define  the  degree  of  elongation,  p  ,  by 

P  =  aT/gx  •  (A- 16) 

Then,  eliminating  and  <j  among  (A-ll),  (A-13)and 

1  (A- 16),  we  obtain 

! 

I 

! 

f  ^  2  7 

— —  =  i  +  l1  ~  P)  (A- 17) 

c  2  2p 

ro 

Finally  we  shall  calculate  the  relative  amount  of  substance 
j  contained  between  r  =  0  and  r  =  n  <?rc  (n  >  0)  for  a  radially 

;  symmetric  distribution  of  a  Gaussian  type. 
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Some  empirical  relations  between  diffusion  characteristics  are 
investigated  by  the  use  of  carefully  examined  data  from  dye  release  experiments 
in  the  surface  layer  of  the  sea.  These  data  cover  a  time  scale  of  diffusion 
ranging  from  1  hour  to  1  month  and  a  length  scale  from  100  m  to  iOO  km.  Two 
"oceanic  diffusion  diagrams"  are  prepared,  one  showing  horizontal  variance 
versus  diffusion  time  and  the  other  showing  apparent  diffusivity  versus  the  scale 
of  diffusion,  The  overall  behaviors  of  the  horizontal  variance  and  of  the  apparent 
diffusivity  are  evidently  different  from  those  which  the  similarity  theory  of 
turbulence  deduces.  However,  there  still  remains  a  possibility  that  the 
similarity  theory  may  be  valid  locally  in  time-  or  length-scale  with  a  variable 
parameter,  i.  e.  ,  the  rate  of  turbulent  energy  transfer.  The  diagrams  provide 
a  practical  means  to  predict  the  rate  of  horizontal  spread  of  aubBtance 
introduced  from  an  instantaneous  point- source  as  well  as  the  apparent  diffusivity 
as  a  function  of  the  size  of  diffusing  patch. 
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